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Fermi inhibition is a quantum statistical analogue for the 
inhibition of spontaneous emission by an excited atom in a 
cavity. This is achieved when the relevant motional states 
are already occupied by a cloud of cold atoms in the internal 
ground state. We exhibit non-trivial effects at finite temper- 
ature and in anisotropic traps, and briefly consider a possible 
experimental realization. 



The inhibition and enhancement of spontaneous emis- 
sion of an atom by embedding it in a cavity has been 
one of the most striking predictions in the field of cavity 
QED Several experiments have and are currently 

being carried out to test this prediction. The explanation 
of this phenomenon is very simple: according to Fermi's 
golden rule, the spontaneous emission rate T is propor- 
tional to the density of modes of the electro-magnetic 
field available at the atomic transition frequency; the cav- 
ity modifies the mode structure, and therefore this rate. 
In particular, if all the normal modes of the cavity are 
far off resonance with respect to the atomic transition 
frequency, there are no modes available where a photon 
can be emitted, and therefore spontaneous emission is 
inhibited. 

In this paper we analyze how the spontaneous emission 
rate T of an atom can be inhibited by using quantum 
statistical principles. The idea is to surround a fermionic 
atom in an internal excited state |e) by a gas formed by 
identical particles in an internal ground state \g). Ac- 
cording to Fermi's golden rule, F depends on the occu- 
pation number of the different motional states to which 
the atom can be transferred after spontaneous emission. 
In particular, if they are already occupied then the atom 
cannot decay. This can be viewed as well as a simple con- 
sequence of Pauli's exclusion principle. Thus, in analogy 
to cavity QED one can say that spontaneous emission can 
be inhibited by reducing the motional "modes" available 
to the excited atom. 

From a qualitative point of view, the inhibition of spon- 
taneous emission can be easily understood. Consider a 
set of N fermions in the state \g) and one fermion in an 
excited state |e) that can only decay into \g). Assume, for 
simplicity, that all the atoms are trapped in a harmonic 
potential at a given temperature. In this problem there 
are three energies scales which determine the properties 
of r. First, ksT, the thermal energy of the particles. 
Second, the Fermi energy Ep of the ground state atoms, 
which is the maximum energy of the fermions at zero 
temperature. Finally, the recoil energy = h^k^/2m 



(where feg is the wavevector corresponding to the tran- 
sition |e) — > \g) and m is the atom mass), which is the 
averaged motional energy acquired by the excited atom 
in the emission process. Typically, the motional states 
with energies E < Ep — ksT are occupied, whereas those 
with E > Ep + ksT are basically available. Thus, if the 
initial motional energy of the excited atom plus the one 
gained in the emission process ~ ksT + Ef; < Ep — ksT, 
spontaneous emission will be inhibited. As the temper- 
ature decreases or the Fermi energy increases (equiva- 
lently, the number of ground state atoms increases) this 
phenomenon becomes more pronounced. A quantitative 
analysis of this phenomenon, on the other hand, is much 
more complicated since in order to evaluate the results 
given by Fermi's golden rule one needs to perform heavy 
numerical calculations. In this paper we have performed 
these numerical calculations, as well as obtained simple 
analytical formulas for some interesting regimes. In par- 
ticular, we have studied the inhibition of spontaneous 
emission of photons along different spatial directions both 
for isotropic and cylindrically symmetric traps as a func- 
tion of the temperature, Fermi energy and anisotropy 
parameter. We have also analyzed a situation which we 
think it is closer to the foreseeable experimental possibil- 
ities. 

We consider an excited atom confined in a harmonic 
potential at some given temperature T. Let us denote by 
ro(r2)c?il the spontaneous emission rate of photons in the 
solid angle dfl around the direction defined by f2 for a sin- 
gle atom, and by r(r2)c?f2 the corresponding value in the 
presence of N ground state fermions confined in the same 
potential at the same temperature T. We are interested 
in the "modification factor" Mf{il) = T{il)/To{^}) which 
gives us information about the modification of sponta- 
neous emission along a given direction f2. Using Fermi's 
golden rule we obtain 
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Here, = (g(7ic^/feBT)(?.n-M) + are the Fermi- 

Dirac and P,n = Ppe^^'^'^/''^^)'^'" the Boltzmann dis- 
tribution for the ground and excited state atoms; \n) = 
\nx)\ny)\nz) denote the eigenstates of the harmonic oscil- 
lator with frequencies (wa;, Wy, w^) = luv along the prin- 
cipal axis 1^. The chemical potential and the number 
of atoms N can be expressed in terms of Ep and lu in 
the usual way. Finally, k{fl) /ko is a unit vector along 
the direction given by fl. The interpretation of this for- 
mula is very simple. The probability density that the 
excited atom is transferred from \rh) to \n) by emitting a 
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photon along the direction Q is |(n|e~''^^^^ ''|TO)p. This 
matrix elements has to be multiplied by the probability 
that the initial state and the final state are occupied [P^n 
and (1 — Ffi)] and summed over. Note that, in order 
to determine Mf, one has to perform six infinite sum- 
mations for each value of Q, which in general cannot be 
evaluated even using a supercomputer. In the follow- 
ing we will consider that Wz/A = lUx = cOy, in which 
case Mf{n) = Mf{9) only depends on the polar angle 9 
and the problem becomes tractable numerically. Further- 
more, we will also derive some analytical approximations 
to this formula. 

To begin we consider an isotropic trap (A = 1). In this 
case, using the spherical symmetry one can show that the 
emission rate becomes independent of the photon direc- 
tion and Eq. (|l|) reduces to 
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^ {l + l)\{n-l\e^^^^\m~l)\\ (2) 
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where Tq and T are the spontaneous emission rates in free 
space and in the presence of quantum statistical effects, 
respectively. Now, the states |n) denote one dimensional 
eigenstates of the harmonic oscillator and P„i and Fn are 
the Boltzmann and Fermi distributions at a given energy 
nhu. These sums can be performed numerically; the re- 
sults are shown in Fig. |l| as a function of temperature 
and Fermi energy. One can see clearly that, as predicted 
by the simple arguments given above, Fermi inhibition is 
strongest at low temperatures and large Fermi energies, 
but does persist at relatively higher temperatures. 
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FIG. 1. Alf as a function of temperature, for a trap with 
Er — 25hui. The solid curves represent Fermi energies rang- 
ing in steps of ten from Ef = lOhui (top) to Ef = SOhcj 
(bottom); the dotted curves are the semi-classical approxi- 
mation. The inset shows Mf at T = as a function of Fermi 
energy. 



At zero temperature the Fermi-Dirac distribution be- 
comes a step function, and the sums in Eq. (|^) may be 
evaluated explicitly. The result is an incomplete Gamma 
function Mf = T[Ep/{fiuj), Ej^/{hoj)], whose integral 
representation expresses the basic competition between 
Fermi energy and recoil energy Q| We have plotted this 
result in the inset of Fig. |l|. As long as the Fermi energy 
is much smaller than the recoil energy, the spontaneous 
emission rate of the atom is not substantially different 
from the free case; this is because the matrix element 
in Eq. (|^) is concentrated around n — m ~ -ii{Eji/huj) 
(partial conservation of momentum due to the finite 
strength of the confining potential). So Fermi inhibition 
at zero temperature is negligible as long as levels around 
Eji/{huj) are still unoccupied. Once the Fermi sphere 
expands beyond these levels, however, the probability of 
emission drops sharply, since to emit an optical photon 
and move into an available motional state then requires 
an improbably large impulse from the trapping potential. 

At high temperatures, the sharp edge of the Fermi- 
Dirac distribution erodes; eventually this means that we 
can analyze Eq. (|^) semi-classically. The matrix element 
for the m — > n transition is concentrated around the 
values satisfying fi.u;(n — m) = ±Eji with width 6E ~ 
2(7iwi?_R)^/^\/2m, where mhuj is typically < 2kBT (oth- 
erwise Pm — 0). On the other hand, the Fermi-Dirac dis- 
tribution changes appreciably only at energy scales ksT, 
i.e, Fn ~ Fn' for nu!\n-n'\ < keT. Thus, USE <C ksT, 
i.e. Er < ksT we can replace F„ F„^±ER/(hu,) in 
Eq. (0). The sum over n then equals one; and at high 
temperature we can also approximate the sums over m 
by integrals, obtaining 

Mf = — = e/3('7^-iJ^) iog(l + e-/3(''^-^-)) (3) 
To 

As can be seen from Fig. |^, this approximation is actu- 
ally quite good even for moderately high temperatures. 
As T ^ cxD the semi-classical result becomes exact, and 
shows that Fermi inhibition disappears in the classical 
limit, when the chance of any particular level being oc- 
cupied becomes infinitesimal. By combining this observa- 
tion with the fact that for Ep < Er Fermi inhibition re- 
quires temperatures sufficient to populate the levels near 
Efi, we can explain why the upper curves in Fig. |l| actu- 
ally show local minima at finite temperatures. 

We now consider the anisotropic case, A 7^ 1. It is 
remarkable that for an atom in free space spontaneous 
emission does not depend at all on the shape of the 
trap [^. The presence of ground state fermions in an 
anisotropic trap changes this fact: the spontaneous emis- 
sion rate does depend on the polar angle 9 between the 
emission direction and the trap z-axis. We emphasize 
that this is a specific consequence of quantum statistics. 

Returning to Eq. (|l]) for A 7^ 1, we can again begin 
by considering its simplified form at zero temperature. 
Because of the remaining cylindrical symmetry of our 
trap, sums over the three components of n reduce trivially 
to sums over Ux and Uz only. We find 



2 



Mf{e) = e- 
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E E ^ ^ , [sin ^]^"- [cos g]^"- , (4) 
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where n = max(0, ^^^^""^ ), and wc have used = 
Efi/{huj), where rj is the so-called Lamb-Dicke parame- 
ter. 
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FIG. 2. M/(6l) in a trap with = 25huj and A = 11, 
for (a) T = 0, for Ep/ihui) ranging from 1 (outermost) to 
45 (innermost) in increments of 11; and (6) Ef = 45?iaj for 
T/{hu!) ranging from 0.5 (innermost) to 3.5 (outermost) in 
increments of 0.5. The radial scale is linear, and the same for 
all curves; the radius of the circular frames is one. 

The results of evaluating numerically Eq. are dis- 
played in Fig. 2. Examining these results reveals, in ad- 
dition to fine structure for which no simple, qualitative 
explanation is apparent, a rather surprising general fea- 
ture. As Ep rises from zero, Fermi inhibition is initially 
strongest in the direction in which the trap is stiffer; but 
at high Ep, it is emission in the soft direction which is 
most inhibited. One might think to explain the greater 
initial inhibition in the stiff direction by reference to the 
lower density of states into which the emitting atom must 
recoil, if it recoils in the stiff direction. This tempt- 
ing explanation is incorrect, however, because the larger 
trap force in the stiff direction allows greater violation 
of momentum conservation, and consequently the range 
of energy levels into which the atom may recoil is also 
greater. It is actually trivial to show that the greater 
width of |(n|e*'^^|m)p exactly compensates for the lower 
density of states, if all final states |n) are available (which 
is why emission by a single atom is independent of the 
trap shape). And this reveals the true explanation for the 
crossover from soft to stiff directions. While | {n\e^'^^ |^ 
is peaked at the same energy Ep regardless of A, the 
wider distribution for the stiff direction is first to be en- 
croached on by the rising Fermi level, which prevents 
recoil into occupied states. The wider distribution is also 
last to be completely submerged as Ep rises past Eji, 
so for large Ep spontaneous emission is preferentially in 
the stiff direction. And the crossover point is once again 
Ep — Er. 

While Fig. 2 presents results for a prolate ('cigar- 
shaped') trap, so that 'the soft direction' is actually 
the equatorial plane, to describe an oblate ('pancake- 



shaped') trap one need merely rotate the circular plots by 
ninety degrees. One piece of curious structure in Mf{9) 
as a function of Ep which is not shown in Fig. 2(a) is 
that if A is an integer, M/(0) and A//(7r) are A-fold de- 
generate as functions of integer Ep/{Tiuj). If A is the in- 
verse of an integer, M/(±7r/2) are similarly degenerate. 
Fig. 2(b) shows that at finite temperatures Mf{6) be- 
comes more isotropic. Indeed, at very high temperature 
complete isotropy is restored, as Pauli exclusion becomes 
insignificant in the classical limit. Therefore the occur- 
rence of the spatial anisotropy of the spectrum can also 
be used signaling the onset of the degenerate regime in 
the cooling process of the gas. 

So far we have assumed that the one excited atom 
had initially been brought into the trap without any dis- 
turbance of the equilibrium population of ground state 
atoms. But let us now consider the experimentally more 
straightforward scenario in which the excited atom has 
been created by applying a laser pulse to a ground state 
cloud. To ensure that no more than one excited atom is 
involved, we suppose a weak pulse. On the other hand, 
we will assume that the pulse is short compared to the 
internal and external dynamics of the atoms, so that we 
can consider it instantaneous. Taking the excited atom 
from the ground state cloud introduces two significant 
changes. Firstly, even at zero temperature there is a sig- 
nificant probability of the excited atom having energy 
above Ep, because it may have been excited from a mo- 
tional state near the Fermi level. And secondly, the exci- 
tation of one atom leaves behind a hole in the Fermi sea, 
into which the atom can always recoil. 

Recalculating the rate of spontaneous emission in an 
isotropic trap, from the initial state with one atom ex- 
cited from the cloud as we have described, we obtain 



Mfin) = i + 
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where Mf{VL) = V{Q.) /[NVL]^tlTQ{n)] is the emission rate 
in the direction O given that an atom has indeed been 
excited. VLr is the Rabi frequency and to the duration 
of the applied laser pulse. The Cj^,{Vl) are the transi- 
tion matrix elements between the states \l') and \l). The 
direction J7 is specified by the polar an azimuthal angles 
{0, (j)) with respect to the propagation direction of the ex- 
citing laser pulse. Obviously, the results for the isotropic 
trap depend only on 9. 

The first term within the large brackets in Eq. (|^) 
represents processes in which an atom is excited by the 
laser pulse, then spontaneously emits a photon and re- 
coils back into the 'hole' from which it came. Because 
the final state of the gas is the same no matter which 
atom undergoes this process, absorption and re-emission 
by all N atoms interfere constructively, and so this term 
is of order N'^; but since the atom recoils back into the 
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state from which it was kicked by the laser pulse, the re- 
emitted photon is overwhelmingly likely to be in the same 
direction as the pulse |^ . The second term in the brackets 
of Eq. (6) describes Fermi-inhibitcd spontaneous emis- 
sion just as we have discussed in the first part of this 
Letter, but with the more complicated initial state for 
the excited atom. It is this term that is of interest, but 
since it is only proportional to it is dominated by 
the first term within a cone around the forward direc- 
tion. Evaluating the two terms for small angles from 
the laser direction, one finds the width of this cone to 
be \Se\ ~ (?7iVi/6)-i ^ l/(fcoL), where L denotes the 
atomic cloud size (see Fig. 3). This is the familiar results 
of diffraction of light by a transparent object of size L. 

We have plotted in Fig. 3 Mf{6) as function of 6. For 
9 > 66, Fermi inhibition is indeed significant, though 
it is nowhere close to the near-total suppression found 
in the cases where the excited atom is initially thermal. 
This is because there is always a significant fraction of 
the ground state atoms within two recoil energies of Ep, 
and these can be excited and spontaneously emit without 
much hindrance. 



sions, which is a good approximation since s-wave scat- 
tering is suppressed for spin polarized Fermions. On the 
other hand, Fermi's golden rule neglects the effects of 
dipole-dipole interactions and reabsorption, and there- 
fore our model assumes that the atomic gas is optically 
thin, meaning that the probability of reabsorbtion of the 
emitted photon Preab = t^(-2-) M is negligible. Some 
straightforward calculations show that this implies the 
condition 0.15E'j/{Efih,u!) <C 1, which for stiff traps will 
still allow significant Fermi inhibition. 

In conclusion, we have analyzed the possibility of ob- 
serving Fermi inhibition of spontaneous emission in an 
atomic gas. We have derived simple analytical formulas 
for the modification factor Mj and compared with nu- 
merical results for both isotropic and cylindrically sym- 
metric traps. We have also studied the situation in which 
an atom is excited by a weak short laser. 
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FIG. 3. M/ at T = for a trap with Er = huj and 
Ef = lOhuj. The solid line corresponds to Eq. [s] and the 
dashed line to the case where forward scattering is suppressed 
by a Raman pulse excitation. In both cases the above lying 
curve is missing the Fermi-inhibition term (the last one in Eq. 



To see Fermi inhibition without simply waiting for rare 
emissions outside the forward cone, one possibility is to 
make the exciting pulse a Raman pulse, so that the recoil 
momentum when the atom emits one photon is different 
from its recoil on absorbing the two Raman photons. By 
this means one can ensure that recoil back into the hole 
the excited atom came from is strongly suppressed by 
momentum conservation, and the order N'^ term in Eq. 
(6) can be essentially eliminated. This effect is shown 
by the dotted curves in Fig. 3. We emphasize, however, 
that the Raman laser should involve a dipolc forbidden 
transition in order for the |e) —^ \g) to be allowed, which 
makes the experiment more challenging. 

Finally, we would like to mention that in the calcula- 
tions presented here we have neglected atom-atom colli- 
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